##Exerccios 4 - Anlises Exploratrias
#BIE572
#Francisco M. Ulloa S.

##4.2 Cervejas
cervejas <-c("chope","lata","garrafa","chope","garrafa", "garrafa","lata","lata","nenhuma","lata","garrafa","garrafa", "garrafa","lata","lata","lata","garrafa","lata","chope","nenhuma", "garrafa","garrafa","garrafa","chope","garrafa","garrafa","chope","garrafa","lata","lata")

#4.2.1 Represente este resultado como um grfico de barras e um dotplot (funo dotchart).

barplot(table(cervejas))
x11()

cervejas.todas<-cervejas
cervejas.todas[cervejas.todas=="nenhuma"]<-0
cervejas.todas[cervejas.todas=="chope"]<-1
cervejas.todas[cervejas.todas=="lata"]<-2
cervejas.todas[cervejas.todas=="garrafa"]<-3
class(cervejas.todas)
cervejas.todas<-as.numeric(cervejas.todas)
class(cervejas.todas)
dotchart(cervejas.todas,groups=factor(cervejas))

#4.2.2 Qual tem maior razo dado/tinta?


##Ambos os grficos mostram informaes claras, porem o dotchart possui uma razo dado/tinta maior, 
#que permite uma melhor visualizao dos poucos dados presentes em "cervejas", e que permitiria 
#um menor gasto de tinta


#############################################################################

#4.3 Caixetais

caixeta <- read.csv("caixeta.csv", as.is=TRUE)

#4.3.1 Construa um histograma do dap dos fustes dos caixetais.

## O dap = o cap/pi
pi
names(caixeta)
table(caixeta$cap)
caixeta$dap<-(caixeta$cap)/pi


hist(caixeta$dap)


#4.3.2 Construa histogramas da altura das rvores para os diferentes caixetais ('local')


library(lattice)
require(lattice)
histogram(~ caixeta$h | caixeta$local )

#4.3.3 H diferenas entre as estruturas (distribuio de tamanhos) dos caixetais?

##Segundo o observado no histograma existem aparentemente arvores de maior tamanho em Jureia em relao aos outros lugares.
##Se poderia tambm supor que as florestas dos caixetais em Retiro esto em estagios mais iniciais que em comparao a Chaus.


##############################################################################

#4.4 Eucaliptos

egrandis<-read.table("egrandis.csv",header=TRUE,sep=";",as.is=TRUE)

#4.4.1 Utilize o grfico boxplot para analisar o DAP de rvores de E. grandis em funo das variveis regio (regiao) e rotao (rotacao).

str(egrandis)
boxplot(egrandis$dap)
egrandis$regiao<-as.factor(egrandis$regiao)
egrandis$rotacao<-as.factor(egrandis$rotacao)
boxplot(egrandis$dap~egrandis$regiao)


#4.4.2 Avalie a normalidade da altura do conjunto total de rvores com um grfico quantil-quantil contra a distribuio normal.

qqnorm(egrandis$ht)
qqline(egrandis$ht, col="red", lwd=3)

# Se observa no grfico quantil-quantil que os dados de altura "ht" no seguem uma distribuio normal.

#####################################################################################

#4.5 Mais Caixetais

caixeta <- read.csv("caixeta.csv", as.is=TRUE)

#4.5.1 Analise a relao dap-altura ('dap' e 'h') em funo do caixetal (local) com a funo plot, mas somente para as rvores 2) de caixeta (Tabebuia cassinoides).


caixeta$local<-as.factor(caixeta$local)
caixeta$especie<-as.factor(caixeta$especie)


pi
caixeta$dap<-2*(caixeta$cap/(2*pi))
str(caixeta)
par(mfrow=c(2,2))
plot(dap~h,data=caixeta,subset=especie=="Tabebuia cassinoides"&local=="chauas")
plot(dap~h,data=caixeta,subset=especie=="Tabebuia cassinoides"&local=="jureia")
plot(dap~h,data=caixeta,subset=especie=="Tabebuia cassinoides"&local=="retiro")
par(mfrow=c(1,1))

#4.5.2Para a mesma relao do item anterior, verifique linearidade com a funo scatter.smooth

par(mfrow=c(2,2))
scatter.smooth(caixeta$dap~caixeta$h,data=caixeta,subset=caixeta$especie=="Tabebuia cassinoides"&caixeta$local=="chauas")
scatter.smooth(caixeta$dap~caixeta$h,data=caixeta,subset=caixeta$especie=="Tabebuia cassinoides"&caixeta$local=="jureia")
scatter.smooth(caixeta$dap~caixeta$h,data=caixeta,subset=caixeta$especie=="Tabebuia cassinoides"&caixeta$local=="retiro")
par(mfrow=c(1,1))

#4.5.3 Utilizando o pacote lattice, analise a relao dap-altura ('dap' e 'h') em funo do caixetal (local), mas somente para as rvores 3) de caixeta (Tabebuia cassinoides).

library(lattice)
require(lattice)


xyplot(dap~h|local,data=caixeta,subset=especie=="Tabebuia cassinoides")
